Dependence of the superconducting transition temperature on the doping level in 

single crystalline diamond films 
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Homoepitaxial diamond layers doped with boron in the 10^" - 10^^ cm range are shown to be 
type II superconductors with sharp transitions (~ 0.2 K) at temperatures increasing from to 2.1 
K with boron contents. The critical concentration for the onset of superconductivity is about 5 — 7 
10'^° cm~'^, close to the metal-insulator transition. The H — T phase diagram has been obtained 
from transport and a.c. susceptibility measurements down to 300 mK. These results bring new 
quantitative constraints on the theoretical models proposed for superconductivity in diamond. 

PACS numbers: 74.25.0p, 74.63.C, 73.61.Cw 



Type II superconductivity has been recently reported 
for heavily boron-doped polycrystaUine diamond pre- 
pared either as bulk or thin film samples 0, provid- 
ing a new interesting system for the study of supercon- 
ductivity in doped semiconductors. Based on the great 
bonding strength of the valence band states and on their 
strong coupling to the carbon lattice phonons, various 
theoretical studies 0, have stressed the similarity be- 
tween diamond and the recently discovered MgB2 system 
which shows a surprizingly high Tc value on the order of 
40 K. Those calculations lead to Tc values in the 0.2 to 
25 K range (depending on the boron content) but ignore 
the boron impuritv band 0, 0. An alternative theo- 
retical approach p stresses out the fact that the boron 
concentration range where superconductivity has been 
observed is close to the Anderson-Mott metal-insulator 
transition and suggests an electron correlation driven ex- 
tended s— wave superconductivity in the boron impurity 
band. 

The open questions of the nature of the metal-insulator 
transition (MIT) in diamond and of its correlation with 
superconductivity are of fundamental interest and pro- 
vide ample motivation for this first investigation of the 
dependence of the superconducting transition tempera- 
ture on the doping boron concentration. In this letter, we 
report on magnetic and transport experiments on a set of 
high quality single crystalline epilayers doped in the rel- 
evant 10^" — 10^^ cm~'^ range. We show that Tc rapidly 



increases above some critical concentration 
reaching 2 K for rt^ = 19 10^*^ cm 



-3 



5 - 7 1020 
(see Table 



cm 
l). 

001-oriented type lb diamond substrates were first ex- 
posed to a pure hydrogen plasma. Methane (4%) was 
subsequently introduced and a 0.5/im-thick buffer layer 
of non-intentionally doped material was deposited at 
820° C by the microwave plasma-assisted decomposition 



TABLE I: Sample characteristics : thickness (t), gas phase 
ratio ((B/C)5as), boron concentration (ns) and critical tem- 
peratures (Tc) for the studied diamond epitaxial films 



Sample 


t (/im) 


{'B/C)gas (ppm) 


ns (10^° cm"^) 


Tc (K) 


1 


3.0 


1615 


3.6 


< 0.05 


2 


3.0 


1730 


9 


0.9 


3 


3.0 


1845 


10 


1.2 


4 


2.0 


2200 


11.5 


1.4 


5 


0.15 


2800 


19 


2.1 



(MPCVD) of the gas mixture at a total pressure of 30 
Torr. Finally, diborane was introduced in the vertical 
silica wall reactor with boron to carbon concentration 
ratios in the gas phase ((i3/C)gas)ranging from 1500 to 
3000 ppm. With a typical growth rate of 0.9^m/hr these 
deposition conditions led to 0.1 to 4/xm-thickj?+— type di- 
amond layers. Secondary Ion Mass Spectroscopy (SIMS) 
depth profiles of ^^B~, ^^C" and ^^B^^C~ ions were 
measured using a Cs+ primary ion beam in a Cameca 
Ims 4/ apparatus. As shown in Fig.l for samples 2, 4 
and 5, these profiles were found to be flat or with a slow 
decrease toward the buffer layer. The boron atomic den- 
sities ub shown both in Fig.l and Table 1 were derived 
from a quantitative comparison to a SIMS profile mea- 
sured under the same conditions in a i?— implanted di- 
amond crystal with a known peak boron concentration 
of 2.4 lO'^" cm^'^. For thin enough samples, the profile 
yielded also the residual doping level in the buffer layer, 
around 10^^ cm"'^. Moreover, the single crystal and epi- 
taxial character of the MPCVD layers was checked by 
high resolution X-ray diffraction, yielding shifted narrow 
lines with smaller linewidths at half maximum than for 
the lb substrate (typically 10 arcsec for the 004 Bragg 
spot of the epilayer, instead of 13 arcsec for the sub- 



2 



^ 10- 
o 



C 

o 

u 

B 
o 

o 

o 
pa 



10" 



10" 



10" 



^ #5 


' ^ #2 


> 




#4 




sputtering 
stop 


\ buffer 


(m = 23) \^_„_^ 



50 100 150 200 250 

Sputtered depth (nm) 




FIG. 1: SIMS profiles for ion mass 23 obtained using Cs^ 
primary ions on samples 2, 4 and 5. In the case of the thicker 
samples, sputtering was interrupted before reaching the buffer 
layer. 



FIG. 2: Dependence of the superconducting transition tem- 
perature Tc on the boron concentration ub : closed squares : 
this work, open circle : from [l].The different Tc values were 
obtained from the onset of the diamagnetic signal (see inset 
for the real (closed symbols) and imaginary (open symbols) 
parts of the magnetic susceptibility. 



strate) . The chemical composition , structural 0, 1^ 
and optical 0, S characteristics of these layers at 
room temperature are well documented and have been 
reviewed recently, together with some preliminary trans- 
port measurements 8]. 

The superconducting temperatures have been deduced 
from ac-susceptibility measurements (xac)- The films 
have been placed on top of miniature coils and Tc has 
been obtained by detecting the change in the self induc- 
tion of the coils induced by the superconducting tran- 
sitions. For fully screening samples, we observed a 4% 
drop of the induction L (~ 1 mH). Small ac-excitation 
fields {oj ~ 99 kHz and hac ~ a few mG) have been ap- 
plied perpendicularly to the films. Four Au/Ti electrodes 
were deposited on top of the film with the highest Tc for 
magneto-transport measurements. A very small current 
(~ InA) corresponding to a current density on the order 
of 10"'^ Acm~^ has been used to avoid fiux flow dissipa- 
tion. A standard lock- in technique at 17 Hz was used to 
measure the temperature dependencies of the sample re- 
sistance at fixed magnetic fields. The measurements were 
performed down to 50 mK (xac) and 300mK (transport) 
and the magnetic field was applied perpendicularly to the 
doped plane of the sample. 

The critical temperatures (see Table 1) reported in 
Fig. 2 have been deduced from the onset of the diamag- 
netic signal (see inset of Fig. 2). The susceptibility has 
been rescaled to —1 in the superconducting state and a 
TV ~ 0.9 demagnetization coefficient has been used (the 
values of Tc do not depend on this choice). This choice 
for TV leads to a dissipation peak on the order of 0.2 — 0.3 
as expected for the non linear regime. This is consistent 
with the observation of a slight increase of the width 



of the transition for increasing hac values. As shown in 
the inset of Fig. 2 (and Fig. 3), the transitions are very 
sharp (with a width ^ 0.2 K) allowing an accurate de- 
termination of Tc and stressing out the high quality of 
our samples. In comparison, the polycrystalline samples 
measured in previous reports presented a much larger re- 
sistivity transition width, 1.7 K in 1] and ^ 2.6 K in 

No transition was observed down to 50 mK for the film 
with Ub ^ 3.6 10^*^ cm^'^. For higher boron concentra- 
tions, Tc increases rapidly with doping above some crit- 
ical concentration ~ 5 — 7 10^° cm^'^ reaching 2.1 K for 
n_B = 19 10^" at.cm^'^. The dependance of Tc with dop- 
ing extrapolates towards the data recently obtained by 
[l|. On the contrary, our Tc value for ub ^ 10^^ cm""^ 
(~ 1 K) is much lower than the one recently reported 
by 2\ (~ 4 — 7 K). Our data suggests that these au- 
thors have largely underestimated the boron concentra- 
tion of their polycrystalline samples. They deduced this 
concentration from Hall measurements which are known 
to give results that deviate significantly from the actual 
boron concentration in p"*"— type diamond 7]. It is also 
worth noticing that we observed superconducting tran- 
sitions with Tc on the order of 1 K for boron contents 
~ 0.5at.% whereas recent calculations of the electron- 
phonon coupling led to much smaller Tc values for these 
doping levels [J^fl @ ■ 

The influence of an external field on the superconduct- 
ing transition is displayed in Fig. 3. As shown in panel 
a, the transition is shifted towards lower temperatures as 
the magnetic field is increased. The transition width re- 
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FIG. 3: (a) : temperateure dependence of the electrical re- 
sistance at indicated magnetic fields for the film with n = 



f9.10^" at.cm" 



The dashed lines correspond to the two 



criteria used for the determination of Hc2 (see Fig. 3) i.e. 
R/Rn = 90% and 10% {Rn being the normal state resis- 
tance), (b) temperature dependence of the real part of the 
magnetic susceptibility at indicated magnetic fields for the 
film with n = 9.10^° at .cm . Hc2 has been deduced from 
the onset of the diamagnetic signal. 



mains relatively small up to 1 T and rapidly increases 
for larger fields. In the absence of thermodynamic mea- 
surements, some care should be taken in order to define 
an accurate Hc2 line. This line has been defined from the 
classical R/Rn = 90% criterion (where Rn is the normal 
state resistance). As shown in Fig. 4, the corresponding 
Hc2{T) line can be well described by the classical WHH 
theory |0. We hence get Hc2{0) ~ 1.4 T correspond- 
ing to a coherence length — v^<i>o/27riJc2(0) ^ 150 A 
for ub = 19 10^° cm~'^ ($o being the flux quantum). 
We have also reported on Fig. 4 the line corresponding 
to R/Rn = 10% which gives an indication for the width 
of the transition, pointing out that the transition curves 
rapidly increase for _ff > 1 T. For sample 2, this line has 
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FIG. 4: H-T phase diagram for the films with ns = 19 10^" 
cm~^ and ub ~ 9 10^" cm~"^. The circles (resp. diamonds) 
have being deduced from temperature sweeps of the electrical 
resistance (see fig.2) for R/R„ = 90% (resp. R/Rn = 10%). 
The shadowed area is an indication of the width of the tran- 
sition. The closed squares were defined from the onset of the 
ac-susceptibility (see Fig.2(b)). The full lines are fits to the 
data using the classical WHH theory. 



been deduced from the shift of the diamagnetic response 
with increasing fields. In this case, the rapid broadening 
of the transition is the hallmark of a small critical cur- 
rent density (Jc) again emphasizing the high quality of 
our films. Indeed, in the non linear regime, the suscepti- 
bility is directly related to Jcd/hac (where c? is a charac- 
teristic length scale on the order of the sample thickness) 
and Jc ~ hac/d ^ 1 Acm^^ for x' ^ —0.4. In this case, 
no saturation was observed down to 200 mK, indicating 
some deviation from the classical WHH behaviour. Note 
that for both samples the slope [dHc2/dT]T^o ~ 1 T/K, 
which is almost 2 times smaller than the value reported 
by Ekimov et al. 1\. 

We obtained an almost temperature-independent nor- 
mal state resistivity on the order of 0.5 mfi.cm for ub — 
19 10^*^ cm~'^ suggesting that our samples arc close to the 
metal-insulator transition. On the basis of the criterion 
first proposed by Mott for metal-non metal transi- 

1 /3 

tions, which in its final form {Nc qh = 0.26, where 
an is the Bohr radius) has been verified in a wide vari- 
ety of condensed media 12], the critical concentration in 
p— type diamond is expected to be around iVc = 1 — 2 10^" 
cm~'^ Limitations to this approach arise from the 

discrepancies in the values proposed in the literature for 
the acceptor Bohr radius a//, in line with other inconsis- 
tencies about the valence band parameters of diamond. 
From an experimental point of view, extrapolations [l^ 



4 



of the Pearson-Bardeen model |3l have led to Nc val- 
ues ranging from 1.5 to 3.0 10^° cm~'^, while different 
transport measurements have prompted other authors to 
propose that iVc ~ 7 10^° cm~^ ^^l- It is worth noticing 
that those values have been deduced from room temper- 
ature measurements. To the best of our knowledge, the 
only low temperature estimate {Nc ~ 40 10^° cm~^) has 
been proposed by Tschepe et al. ^ITJ on the basis of 
a scaling analysis of the electrical conductivity. How- 
ever, the absolute conductivity values of their implanted 
samples are much lower than ours despite larger boron 
concentrations . 

The proximity of a MIT in this system, makes it 
very favorable for the observation of quantum fluctua- 
tion effects. The strength of these fluctuations can be 
quantified through the quantity Qu — Reff/Rg where 
Rq — h/e^ ^ 4.1 kil is the quantum resistance and 
Ref / = Pn/s (s is a relevant length scale for these fluc- 
tuations nil). Taking pjv ~ 5 lO^'^fi.cm and s ^ ^(0) ^ 
150A, we obtain a large Qu ratio '^0.1 indicating that 
quantum effects may be important in this system. These 



quantum fluctuations may give rise to the melting of the 
flux line lattice and thus lead to the rapid broadening 
of the resitive transitions observed at low T and large 
H. Another indication for such quantum effects, is the 
almost temperature-independent mixed state resistivity 
above 1.2 T as previously observed in other systems with 
similar (3„values p!g.l2g| 

To conclude, we were able to prepare highly homoge- 
neous and well characterized boron-doped diamond films 
in the 10^" - 10^^ cm~'^ range where superconductivity 
occurs. The value of the critical concentration for the 
onset of superconductivity is on the order of 5 — 7 10^° 
cm~'^. Boron-doped diamond is an ideal system to study 
the occurence of superconductivity close to the metal- 
insulator transition. As a consequence, quantum effects 
are expected to play a significant role as suggested by the 
large quantum resistance ratio Qu ~ 0.1. 

We are grateful to F. Pruvost who grew some of the 
diamond layers and to Dr L. Ortega (Cristallographie, 
CNRS, Grenoble) for the X-ray diffraction experiments. 
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